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ABSTRACT 

We obtain lower limits on the amplitude of convective velocities in the deep solar convection zone 
based only on the observed properties of the differential rotation and meridional circulation together 
with simple and robust dynamical balances obtained from the fundamental MHD equations. The 
linchpin of the approach is the concept of gyroscopic pumping whereby the meridional circulation 
across isosurfaces of specific angular momentum is linked to the angular momentum transport by the 
convective Reynolds stress. We find that the amplitude of the convective velocity must be at least 
30 m s _1 in the upper CZ (r ~ 0.95-R) and at least 8 m s" 1 in the lower CZ (r ~ 0.75i?) in order 
to be consistent with the observed mean flows. Using the base of the near-surface shear layer as a 
probe of the rotational influence, we are further able to show that the characteristic length scale of 
deep convective motions must be no smaller than 5.5-30 Mm. These results are compatible with 
convection models but suggest that the efficiency of the turbulent transport assumed in advection- 
dominated flux-transport dynamo models is generally not consistent with the mean flows they employ. 



1. INTRODUCTION 

Our modern understanding of solar internal dynamics 
rests heavily on global and local helioseismology, with 
perspective provided by theoretical and numerical mod- 
els. This perspective is particularly important in the 
deepest regions of the convection zone where mean flows 
are established and where helioseismic probing is most 
challenging. 

The most reliable result from helioseismology with re- 
gard to the dynamics of the deep convection zone (CZ) 
continues to be the solar inter nal rotation profile ob - 
tained from global inversions ()Thompson et all [2003T) . 
Such inversions indicate that the monotonic surface dif- 
ferential rotation (~ 30% decreases in angular velocity fi 
from equator to pole) persists throughout the convection 
zone with little radial variation. Further dynamical clues 
come from estimates of the meridional flow near the sur- 
face inferred from local helioseismic inversions and pho- 
tospheric observations ((Gonzalez- Hernandez et all 120061: 
Basu fc Antidl2010trUlrichll2010l:lHathawav fc Rightmirel 
2010UHathawavl 1201 1ft . Although various measurement 
techniques can yield disparate results for the radial and 
temporal dependence of the meridional flow, all tech- 
niques generally agree that there is a persistent poleward 
flow near the surface (r > 0.95i?, where R is the solar ra- 
dius) from the equator up to latitudes of at least 60-70°, 
with an amplitude of 10-20 m s _1 . 

Less is known about the structure or amplitude of con- 
vective flows in the deep convection zone. Photospheric 
observations are dominated by solar granulation, with a 
velocity scale of ~ 2 km s _1 and a size scale of ~ IMm, 
but deep convective motions are thought to be much 
slower and larger (e.g. iMiesch fc Toomrej[2009l) . Closer 
scrutiny of the photospheric velocity power spectrum ob- 
tained from Doppler measurements shows a prominent 
peak near a spherical harmonic degree £ ~ 120, cor- 
responding to a size scale of 30-35 Mm and a spectral 



velocity amplitude of 7-8 m s 1 (Hathaway et al 



spectral 
III). 



However, this observed spectral amplitude is significantly 
lower than the true amplitude of the convection due 
mainly to projection effects; the Doppler measurements 
trace mainly horizontal flows near the limb. Taking pro- 
jection and other sampling effects into account (e.g. spa- 
tial and temporal filtering), Hathaway et al. were able to 
match the observed spectrum near £ ~ 120 with a super- 
granular flow component having a typical velocity ampli- 
tude of about 300-400 m s _1 , in agre ement with the cor- 
relation tracking of surfac e features ([DeRosa & Toomrel 
[2001 IRoudier et alll20ll . The power declines steadily 
toward lower £, where one would expect to see signatures 
of deep convection (giant cells). Thus, either giant cells 
have a lower amplitude than supergranulation, or they 
d o not imprint th r ough to the photosphere (or both). 

Hanasog e et al.l (|2010t I2012D have recently searched for 
signatures of subsurface convection in local helioseismic 
inversions at depths of 0.92R, 0.95R, and 0.96iJ. Their 
analysis indicates that the spectral velocity amplitude 
at these depths is less than 10 m s _1 and possibly less 
than 1 m s^ 1 for convective motions with £ < 60 and 
correlation times longer than 96 hours. This is one to two 
orders of magnitude lower than suggested by convection 
simulations and mixing length theory (see £ )4.2I) . 

In this paper we argue that these results, if confirmed, 
would have important implications for the maintnance 
of mean flows in the Sun. In particular, we use funda- 
mental physical arguments together with solar observa- 
tions to derive lower bounds on the amplitude of convec- 
tive motions in the deep solar interior. These estimates 
are based on robust dynamical balances deduced from 
the equations of magnetohydrodynamics (MHD). They 
are not based on mixing length theory and they do not 
depend on any results from numerical simulations, al- 
though they are consistent with both ( ^4.21) . We find 
that the amplitude of convective motions in the upper 
convection zone (r ~ 0.95) must be at least 30 m s^ 1 in 
order to be capable of sustaining the mean flows inferred 
from helioseismology. Similar estimates for the lower CZ 
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(r ~ 0.75-R) imply that convective motions there must 
be at least 8 m s _1 . Reco nciling these lower limits with 
the upper limits found by IHanasoge et al.l (|2010t 12012ft 
is a challenge for global convection models but may be 
possible if the motions responsible for maintaining mean 
flows span multiple scales, with significant power above 
£-60. 

In fJ21 we discuss how differential rotation and merid- 
ional circulation are maintained in the Sun by means of 
the convective Reynolds stress, baroclinicity, and their 
own inertia. This provides a link between convection and 
mean flows that we exploit in ^3] to obtain lower limits on 
the amplitude of convective motions in the deep convec- 
tion zone. In 2] we address the length scale of convective 
motions and the implications of our velocity and length 
scale estimates for helioseimic probing, flux-transport dy- 
namo models, and convection models. Section[S]is a sum- 
mary of our principle results and conclusions. 

2. CONVECTIVE ORIGINS OF MEAN FLOWS 

2.1. Gyroscopic Pumping 

The conservation of angular momentum for a statis- 
tically steady flow i n a rotating spherical she ll may be 
expressed as follows dMiesch fc H indman 2011, hereafter 
MH11): 

( P v m )-VC = F . (1) 

where p is the mass density, C = A 2 S1 is the specific 
angular momentum, fi is the rotation rate, A = rsin# 
is the cylindrical radius, and (v m ) = (v r ) r + (vg) 6 is 
the mean meridional flow. This equation is derived by 
averaging the zonal momentum equation over longitude 
and time, with averages denoted by angular brackets (). 
Spherical polar coordinates (r, 9, <fr) are used throughout. 

The J- term in eq. ([T]) is a net axial torque that includes 
contributions from the Reynolds stress, the Lorentz force, 
and the viscous diffusion. Explicit expressions are given 
in MH11. The molecular viscosity is small in stars so 
viscous diffusion can be safely neglected. Neglecting the 
Lorentz force is less justified. However, the magnetic 
pressure gradient averages out when taking the mean so 
the only contribution of the Lorentz force to T is from 
magnetic tension. This generally te nds to suppress ro- 
tational shear (e.g. IBrun et al J 12004) so it is unlikely to 
be the dominant factor in establishing the solar differen- 
tial rotation. In this paper we wish to establish a lower 
limit on how strong the Reynolds stress must be in order 
to balance the advection of angular momentum by the 
meridional flow; that is, the term on the left-hand-side 
(LHS) of ([1]). If it must also balance magnetic tension, 
then the Reynolds stress must be even stronger. There- 
fore, the inclusion of the Lorentz force is not likely to 
change the lower bounds we establish in 

Thus, if we assume the Reynolds stress is the dominant 
contribution to J-, Eq. {jXJ) becomes 

(pv m ) -V£ = -V- ((pAv^» . (2) 

where primes indicate fluctuations about the mean. 
Since V£ is cylindrically outward in the Sun, the im- 
plication here is that the angular momentum transport 
by the Reynolds stress (RHS) establishes the meridional 
flow by inducing a flow toward the rotation axis when 
the RHS is negative (divergence) and away from the 



rotation axis when the RHS is positive (convergence). 
This mechanism is known as gyroscopic pumping (MH11; 
lMcIntvrelll998ft and is suppor ted by convection simula- 
tions and mean- field mo dels (jRempell 120051 : IBrun et"al"1 
[20TTt iFeatherstone et alJl2012f >. 

The fundamental p hysics behind gyros copic pumping 
is disc ussed by MHll. lHavnes et al.l (jl99lh and lMcIntvrel 
(1998). In short, an axial variation in the net torque 
dT ' jdz establishes an axial shear dil/dz which in turn in- 
duces a meridional flow through the Coriolis force f £12.21) . 
In order to determine the nature of the differential rota- 
tion that will ultimately be established, one must con- 
sider the meridional force balance as we do in the next 
section. However, the equilibrium structure of the merid- 
ional flow is surprisingly less sensitive to the meridional 
force balance and is instead regulated mainly by eq. ([T]). 
This is because the primary contribution to V£ is from 
the mean (globally averaged) rotation rate of the Sun, 
£7o (see MH11 and §3.21 below), so changes in the differ- 
ential rotation, Vf2, do not change the balance in eq. ([IJ 
appreciably. 

2.2. Meridional Force Balance and the Role of 
Baroclinicity 

In ij2.ll we argued that the mean meridional circula- 
tion in the solar envelope is established and maintained 
through gyroscopic pumping, as expressed by eq. ((XJ) . 
Since the process of gyroscopic pumping is mediated by 
the Coriolis force, this is equivalent to saying that the 
meridional circulation is maintained by the inertia of the 
differential rotation (MH11). However, it is well known 
that thermal gradients can also establish meridional cir- 
culation by means of baroclinicity and it is sometimes 
argued that the meridional flow in the solar convection 
zone may be baroclinic in nature. In this section we ar- 
gue that this is not the case. In particular, we consider 
the meridional force balance in the solar convection zone 
within the context of solar observations and we argue 
that gyroscopic pumping, as expressed in eq. (JI]), rather 
than baroclinicity, is the primary mechanism by which 
the solar meridional circulation is maintained. 

The ambiguity arises because the thermal energy equa- 
tion can be expressed in a manner analogous to eq. ([TJ: 

(pv m ) V (S) = Q . (3) 

where Q involves the (negative) divergence of the convec- 
tive entropy flux, the radiative diffusion, and the viscous 
and ohmic heating, although the latter are negligible in 
stellar interiors. 

It is clear from equations ([1} and ^ that both me- 
chanical and thermal forcing, T and Q, can induce a 
meridional circulation. Indeed, this is a classical prob- 
lem in the theory of planetary and stellar atmospheres 
(jEliassenl [19511 iReadl ITiM iTassoull Il978l ). In a solar 
context, eq. ([3]) plays an importa nt role in the radiativ e 
spreading of the solar tachocline ()Spiegel fc Zahnl ll992'). 

We can be reasonably certain from solar observations 
that both T and Q are nonzero. If this were not the case, 
then a steady state would only be possible if either C or 
(S) were constant on streamlines of the meridional mass 
flux. Helioseismic rotational inversions together with ob- 
servations of poleward mass flux in the solar surface lay- 
ers rule out the former (£ constant on streamlines). The 
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latter is ruled out by mass conservation ((ty) must be 
nonzero in the deep convection zone to sustain the pole- 
ward mass flux in the surface layers) and the requirement 
that at least some portion of the convection zone be su- 
peradiabatic d (S) jdr < 0. 

The relative contribution of mechanical and thermal 
forcing to establishing the meridional flow becomes more 
clear if we consid e r the mean zonal vorticity equation 
(e.g. iMieschl 120051: iBalbus et al.l 120091 MH11). This is 
obtained from the meridional components of the MHD 
momentum equation, averaged over longitude and time. 
Previous work based on numerical convection simulations 
(iBrun fc Toomrel l200l iMiesch et all l200ll. mean-field 
models jKitchatinov fc: RudigeHll995tlRempelll2005f) . and 
theoretic al interpretation of helioseismic rotational in- 
versions ([Balbus et al.l [20091 ) suggest that the dominant 
contributions to this equation are the Coriolis and cen- 
trifugal terms associated with the differential rotation 
as well as the zonal component of the baroclinic vector. 
The latter is proportional to the mean latitudinal en- 
tropy gradient if the stratification is nearly hydrostatic 
and adiabatic and if the equation of state is that of an 
ideal gas. Thus, neglecting the meridional components of 
the Reynolds stress, the Lorentz force, the viscous diffu- 
sion, thermal fluctuations, and quadratic terms in (v TO ) 
yields 



d_ 

di 



A 



on 2 

dz 



9 d(S) 
rC P 88 



(4) 



where u)<p is the longitudinal component of the fluid vor- 
ticity ijj = Vxv, S is the specific entropy, g is the gravi- 
tational acceleration, and Cp is the specific heat at con- 
stant pressure. 

We have retained the time dependence in eq. Q in 
order to illustrate how mechanical and thermal forcing 
influence the evolution of the meridional flow. If the 
balance in equation ([1]) is not satisfied, this will lead to 
a change in the specific angular momentum dC/dt = 
X 2 8il/8t that will in turn influence the meridional flow 
through the first term on the RHS of eq. Likewise, 
any imbalance in equation @ will change the baroclinic 
(second) term through 8 (S) /8t. 

The system will evolve nonlinearly toward thermal 
wind balance (TWB), described by an equation that is 
now well known: 



dtl 2 

dz 



9 d(S) 
rXCp 88 



(5) 



Most current theoretical and numerical models of so- 
lar mean flows attribute the conical nature of the so- 
lar differential rotation profile (9S1 2 /dz ^ 0) to latitu- 
dinal entropy gradi ents through eq. ([5]). In particular, 
IBalbus eT al. (2009) have argued that the orientation of 
the isosurfaces in the solar CZ inferred from helio- 
seismology follows from thermal wind balance, Eq. ([5]), 
together with the hypothesis that isorotation and isen- 
tropic surfaces coincide. Furthermore, the existence of 
the near-surface shear layer suggests that there is a tran- 
sition near r ~ 0.95-R below which the Rossby num- 
ber is less than unity and thermal wind balance prevails 
(MH11). 

Although we agree that it plays an essential role in 
determining the orientation of the angular velocity iso- 
surfaces in the solar convection zone, we note that baro- 



clinicity due to axisymmetric thermal gradients cannot 
account for the existence of the solar differential rota- 
tion. The zonal component of the Reynolds stress must 
also contribute. This is demonstrated in the Appendix. 

Here we wish to focus on the approach to thermal 
wind balance by means of the time-dependent vorticity 
equation (j4]). Although the meridional circulation itself, 
represented by (w^), drops out of the balance equation 
in a steady state, any imbalance between j the inertial 
and baroclinic terms on the RHS will induce a flow. 
Consider for example the northern hemisphere (NH). 
There it is we ll known from helioseim ic observations that 
8Vt 2 /8z < ([Thompson et al. 2003). This will tend to 
induce a counter-clockwise (CCW) circulation ((w^) < 0, 
poleward in the upper convection zone) that will act to 
make the f2 profile more cylindrical (dtt/dz — 0) in ac- 
cordance with the Taylor-Proudman theorem. In order 
to oppose this and achieve TWB [Eq. J§J], there must be 
a poleward entropy gradient (d (S) /dd < in the NH). 
In short, the inertial term tends to induce a CCW cir- 
culation in the northern hemisphere while the baroclinic 
term tends to induce a clockwise (CW) circulation (vice 
versa in the southern hemisphere). 

In the approach to equilibrium, one of these terms must 
act as the driver, accelerating the meridional flow, while 
the other acts as a resistance, opposing the acceleration 
until a balance is achieved. The observed poleward sense 
of the meridional flow near the solar surface (SJTJ) indicates 
that the driver is in fact the inertial term cx 80,/dz. As 
demonstrated in £|2. 11 this links the mean flows directly to 
the convective Reynolds stress and enables us to estimate 
the amplitude of the convective velocities based on the 
observed differential rotation and meridional flow. 

There is one po tential caveat to this conclusion. The 
recent analysis bv lHathawavl (|2011[ ) based on autocorre- 
lation of photospheric Dopplergrams suggests that there 
may be a reversal of the meridional flow at r ~ 0.95i?. 
If this is indeed true, than we cannot rule out a CW 
circulation cell (NH, CCW in the south) in the deep 
convection zone driven by baroclinic forcing (implying 
poleward flow near the base of the CZ). However, this 
is in conflict with helioseismic inversions which suggest 
the meridional flow remains poleward well below 0.95-fi 
(IGiles et alJll997t IBraun fc Fanlll99l IChou fc Daill2001fc 
iBeck et al.ll2002D . Furthermore, poleward flow near the 
base of the CZ would be disasterous for flux-transport 
dynamo models ( 34.3[) . This in itself does not preclude 
the presence of such a cell but it demonstrates that pole- 
ward flow at the base of the convection zone is contrary 
to our current understanding of solar interior dynamics. 

3. AMPLITUDE OF CONVECTIVE VELOCITIES 

3.1. Fundamental Expressions 

In this section we use Eq. @ to obtain lower limits on 
the amplitude of convective velocities throughout the so- 
lar convection zone. The motions we are most interested 
in are the motions that are responsible for establishing 
the solar differential rotation by means of the Reynolds 
stress. Thus, unlike solar granulation, they must be large 
enough and slow enough to sense the rotational influence 
and spherical geometry. The mere existence of the solar 
differential rotation is evidence for the presence of such 
motions. In what follows, we will quantify what we mean 
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by large enough and slow enough. 

We begin by estimating the amplitude of the Reynolds 
stress inside the divergence operator on the RHS of Eq. 



(pAv>;) ~ pXeV c 2 



(6) 



where V c is the characteristic amplitude of the convec- 



tive velocity and e = 



V r . 2 is a correlation coef- 



ficient describing the efficiency of the convective angular 
momentum transport. If flows are perfectly correlated, 
the transport is efficient and e = 1. However, in practice 
e will fall somewhere between zero and one. Convection 
simulations such as that shown in Fig. [4] below yield e ~ 
0.1-0.2. Similar values of are g i ven b y the mean-field 
theory of IKitchatinov fc Riidigerj (|2005D . There the am- 
plitude of the non-diffusive component of the Reynolds 
stress, A, is of order vt^oH where vt ~ V c 2 t c /3 is the 
turbulent viscosity, r c is the correlation time, and T-L is 
a nondimensional normalization factor that depends on 
the inverse Rossby number Vt* — 2r c 51o- This yields 
e ~ AV^" 2 ~ il*T~L/6, which according to their model, 
is about 0.1-0.2 in the lower CZ. 

Note that we use a single velocity scale V c to charac- 
terize the convection but this does not necessarily im- 
ply that the velocity field is isotropic. On the contrary, 
the velocity scale is defined in terms of the Reynolds 
stress through equation (|6]), which requires some degree 
of anisotropy. Thus, V c can be regarded in general as 
the geometric mean of the two velocity components that 
dominate the Reynolds stress. 

Throughout the paper we assume that the density p 
is sphericall y symmetric and is approximately given by 
Model S of IChristensen-Dalsgaard et all (|1996t ). Then 
substituting ((6]) into ([2]) and noting that < e < 1 pro- 
vides a lower limit on V c 



\: - [ - v m \vc 



1/2 



> 



(5V m \VC\) 



1/2 



(7) 



where V m is a characteristic amplitude of the meridional 
flow and S — L t /\ where L t is the length scale associated 
with the divergence operator in ([2]) and A = r sin 9 as 
above. 

As above with the convection, we emphasize that the 
meridional flow is not isotropic. This is particularly 
true in the upper and lower boundary layers of the CZ 
where the meridional flow is predominantly horizontal 
and at the equator where it is predominantly vertical. 
The proper interpretation for V m is the amplitude of the 
flow across £ isosurfaces, which are approximately cylin- 
drical. So, V m should be regarded as the flow toward 
or away from the rotation axis, which may be predomi- 
nantly radial or latitudinal depending on the latitude. 

A useful variation of equation ([7]) can be obtained if 
we consider only the component of £ involving the mean 
rotation rate fio- This gives |V£| ~ 2Vh where Vn — 
Afio is the zonal velocity associated with the rotation of 
the star. Substituting this into Eq. (0 gives 



V e > (2SV m V n ) 1/2 



(8) 



Thus, if (25) 1//2 is of order unity, then a lower limit on 
the convective velocity can be obtained by taking the 
geometric mean between the meridional flow speed and 



the rotational speed, relative to an inertial frame. Taking 
nominal values of V m ~ 2-10 m s _1 , Qq = 2.7 x 10 
s - \ and A ~ 0.85 R/y/2 would then yield V c > 47-106 
m s . However, as we will see below, (2S) 1 / 2 is likely 
to be somewhat less than unity, bringing these estimates 
down. Still, the extension from (J7J to (JSJ is robust, since 
the differential rotation serves to steepen the £ gradient 
relative to the uniform rotation value of 2Vh (see ^3 . 2[) . 

It now remains to estimate S. In this context, it is 
important to emphasize that L t reflects the scale of tur- 
bulent transport, which is not necessarily the same as the 
scale L c of the convective motions themselves (the latter 
is addressed in £|4.1[) . If the balance expressed by Eq. 
([2]) is to be realized, then this scale must be intimately 
linked to the structure of the meridional circulation. In- 
deed, if our picture of gyroscopic pumping is correct then 
the structure of the meridional circulation is largely de- 
termined by the spatial variation of the Reynolds stress, 
as reflected by L t . 

To illustrate this relationship, consider the simplest 
case of a single meridional circulation cell in each hemi- 
sphere, with poleward flow in the upper CZ and equa- 
torward flow in the lower CZ. Then the LHS of Eq. © 
would be correspondingly negative and positive in the 
upper and lower CZ respectively. To balance this, the 
RHS must also change sign, suggesting L t ~ D/2, where 
D is the depth of the convection zone. If T were to 
exhibit multiple sign changes across the CZ, this would 
produce multiple layered circulation cells in radius such 
that L t ~ D/(2N C ) where N c is the number of cells. For 
the moment arm we can consider mid-latitudes in the 
mid-convection zone, so A ~ r m /\2, where r m ~ 0.85R. 
Thus, we have 



S = 



V2L t 



D 



0.25 



V2r m N c N c 



{V m = | (v e ) | at 45° lat) 
(9) 

where we have used the value for the Sun of D/R ~ 0.3. 
As noted, this corresponds to poleward or equatorward 
flow at mid latitudes in one or more circulation cells. So, 
when used with Eq. ([7]) or © the value of V m should 
correspond to the latitudinal flow speed. 

The caveat that can potentially limit the applicability 
of eq. @ is that, although is regulated by the number 
of nodes in T (sign changes in the Reynolds stress diver- 
gence), the flow amplitude V m is regulated by local gra- 
dients. This is particularly important at the base of the 
CZ where a strong convergence of the angular momen- 
tum flux over a narrow layer could drive an arbitrarily 
strong equatorward flow V m for a given V c . 

To demonstrate this concept, consider a single-celled 
meridional circulation profile (N c = 1) as illustrated in 
Figure [T] We define the radius ro as the turnaround 
radius, where the poleward flow in the upper CZ transi- 
tions to the equatorward flow in the lower CZ. We also 
postulate the existence of a radius r\, below which the 
amplitude of the meridional circulation becomes negli- 
gible. The existence of such a radius can be inferred 
based on the presence of the solar tachocline as deduced 
from helioseismology. In the absence of other forces, the 
combined effects of gyroscopic pumping, thermal wind 
balance, and radiative diffusion would induce a merid- 
ional circulation below the convection zone that would 
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Fig. 1. — Schematic diagram illustrating ro and rj, for an idealized 
single-cell meridional circulation profile. Both radii are indicated 
by dashed lines, although the latter is difficult to distinguish from 
the streamlines of the mass flux, indicated by solid lines. The 
northern hemisphere is shown and the sense of the circulation is 
counter-clockwise. The turnaround radius ro marks a change in 
sign of the latitudinal flow while r b marks the radius below which 
the circulation is negligible. The radius of the Sun is denoted by 
R. Thus, equatorward flow occurs for rj < r < ro and poleward 
flow for ro < r < R. 

wipe out the radial shear in the tachocline, causing the 
latitudinal shear to s pread downward on a time scale of 
several billion years (jSpiegel fc Zahn 1992j). Helioseimic 
inversions indicate that this has not happened, suggest- 
ing the presence of some unknown tor ques that act to 
confine the tachocline (e.g.[Micsch 2005). Thus, the thin- 
ness of the tachocline inferred from helioseismic rotation 
inversions suggests that the amplitude of the meridional 
circulation likely drops off rapidly below th e base of the 
convection zone (|Garaud fc Arreguinll2009l ). 

So, in light of eq. (p}, r separates the the region of 
poleward flow and flux divergence (F < 0), ro < r < R, 
from the region of equatorward flow and flux conver- 
gence (J 7 > 0), ^ < r < r . Equation (J9]) effectively 
assumes that the width of these layers is comparable, 
such that ro — r?, ~ i? — ro. In other words, this equation 
breaks down if the equatorward flow is confined to a much 
smaller region than the poleward flow, or vice versa. In 
this case the effective L t = r$ — r& could become arbi- 
trarily small, implying small values of 5 (<C 1). A similar 
caveat also holds for the multi-celled case N c > 1. 

The mechanical forcing required to establish an asym- 
metric meridional circulation profile such as that de- 
picted in Figure [T] via gyroscopic pumping may or may 
not be convective in origin (asymmetric in the sense that 
ro — n, << R — ro). A localized convergence of angular 
momentum flux could arise from the convective Reynolds 
stress as convective plumes are rapidly decelerated by the 
subadiabatic stratification near the base of the CZ. Alter- 
natively, it could arise from whatever non-convective pro- 
cesses may be involved in tachocline confinement, such 
as large-scale Lorentz forces or Reynol ds and Maxwell 
stresses induced by MHD instabilities (iSpiegel fc Zahnl 
[1991 IGough fc Mclntvrel fl998L IMieschl I2005T I . What- 
ever its origin, the confinement mechanism would have 
to transport angular momentum poleward in order to 
prevent the radiative spreading of the tachocline. This 
would imply a convergence of angular momentum at mid 



to high latitudes that would in turn induce a progradc 
torque and an equatorward flow by means of gyroscopic 
pumping. This may plausibly be confined to a thin 
boundary layer in the vicinity of the tachocline. 

Regardless of how thin the region of equatorward flow 
may be, conservation of mass requires that there be a 
net poleward flow in the remainder of the CZ, implying 
a flux divergence. Might this also occur in a thin layer, 
such that 5 << 1? If so, wouldn't equation ([7]) imply 
that V c could be arbitrarily small? In other words, could 
weak convection maintain the observed mean flows in the 
solar interior by establishing a Reynolds stress that is 
nearly divergenceless in the bulk of the CZ, with narrow 
regions of divergence and convergence in the upper and 
lower boundary layers? In the remainder of this section 
we argue that the answer to this question is likely to be 
no. Although this scenario is in principle consistent with 
the physics of gyroscopic pumping, it can be largely ruled 
out by helioseismic measurements. 

We begin with a direct estimate of 6 in the surface 
layers obtained from helioseismic determinations of the 
meridional flow. As noted in <JIJ local heloiseismic in- 
versions suggest that the radius at which the poleward 
surface flow reverses sign is no shallower than 0.95-R, so 
R — ro > 0.05.R. If we equate this with L t , this provides 
a lower bound for S in the surface layers: 

6 ~ V2L t /R > 0.07 (surface layers) . (10) 

Thus, the factor of 5 1 / 2 in eq. (JT]) is likely to be no smaller 
that 0.26 in the surface layers. 

We now consider the nonlocal nature of gyroscopic 
pumping, which provides a link between regions of 
Reynolds stress convergence and divergence even if they 
are spatially separated (MH11). Local forcing can in 
principle induce a global meridional flow that couples 
the two regions, even if they are localized in the upper 
and lower boundary layers. Furthermore, since isosur- 
faces of C are approximately radial at high latitudes, a 
single circulation cell with poleward flow at the surface, 
downward flow near the poles, and equatorward flow at 
the base of the CZ could in principle be sustained even if 
the high-latitude Reynolds stress were divergenceless in 
the bulk of the CZ {F « 0). 

However, T cannot be zero in the bulk of the CZ at low 
latitudes. Mass conservation requires a radially outward 
meridional flow that must be balanced by the Reynolds 
stress. If it were not, the circulation would quickly ho- 
mogenize C at the equator, establishing an inward Q gra- 
dient (fi oc r -2 ) on a crossing time scale ~ D/V m (about 
6 years for V m ~ 1 m s -1 ). Thus, the £1 profile and 
the horizontal divergence of the meridional flow inferred 
from solar observations provide a robust diagnostic of 
the Reynolds stress at the equator. We now exploit this 
diagnostic to obtain alternative estimates for S and V c . 

Symmetry requires that the radial component of the 
flux divergence on the RHS of {5| vanish at the equator, 
so L t must correspond to the latitudinal convergence of 
the angular momentum flux that sustains the radially 
outward V m . Since the convergence of T must span all 
latitudes where the flow is outward, we can estimate L t 
from the outward mass flux M — pV m 2irrLt. We have 
defined M to be the outward mass flux in one hemisphere 
because in £|3.2l we will equate it to the observed poleward 
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Fig. 2. — Shown is the magnitude of the specific angular momen- 
tum gradient | V£| inferred from helioscismic rotational inversions. 
For the corresponding profiles of Q and C see Fig. 1 of MH11. These 
results are based on RLS inversions of GONG data fro m four non- 
overlapping intervals in 1996, provided by R. Howe IIHowe et al.l 
[20001 ; ISchou et al.ll2002T) . The dotted line indicates the base of the 
convection zone. 

mass flux in the surface layers. Setting A = r at the 
equator then gives us an alternate estimate for 8 

^ = T= t/o 2 <y m = | (v r ) | at 0° lat) . (11) 
Substituting this expression into ([7]) yields 

/ ' \ 1//2 

Vc > [ VZ I (^n = I («r) I at 0° lat). (12) 

Note that the meridional flow only appears in Eq. (|12[) 
through the mass flux JVf, which we can estimate from he- 
lioseismic inversions ( §3.2p . Since C and p are also known 
from helioseismology and solar structure models, Eq. (|12[) 
provides a robust lower limit on the convective velocity 
near the equator that is independent of 5. Given the 
obser ved uniformity of the solar irradiance (Ra st et al.l 
2008), velocity amplitudes at higher latitudes are likely 
to be comparable. 

In summary, Eq. (|12l) provides a reliable lower limit on 
the convective velocity amplitude at low latitudes based 
firmly in observable quantities. The estimate is most reli- 
able at a radius of r ~ 0.95i? (above which we know that 
the meridional flow is poleward) but it can be extended to 
lower radii if a single-celled profile (N c = 1) is assumed. 
Equation (JTUJ) is also robust, depending only on helioseis- 
mic inversions in the surface layers. Together with Eqs. 
10 and (0 it provides a reliable, though somewhat con- 
servative, limit on the convective velocities above 0.95i?. 

Equation (0 can be used together with Eqs. (0 and 
(0 to provide estimates of V c in the lower CZ at mid- 
latitudes. However, these estimates break down if there 
is a large difference in the filling factor between polar and 
equatorial flows. For example, if N c — 1 and if the return 
equatorward flow near the base of the CZ is confined to 
a very small region such that tq — <C R — tq , then the 
effective 5 near the base of the CZ could be significantly 
smaller than the estimate given in Eq. (0. This would 
in turn reduce the lower limit on V c near the base of the 



CZ according to Eqs. (0 and (0. 

In the next section ( i)3.2[) we quantify our estimates by 
turning to photospheric observations and helioseismic in- 
versions. However, before proceeding, we briefly address 
the possibility of multiple circulation cells in latitude. 

In short, the analysis and conclusions of this paper are 
insensitive to the high-latitude structure of the merid- 
ional flow. Our focus is not on the mean flows them- 
selves, but on using mean flow measurements as tools 
to probe the convection. Since helioseismic inversions of 
differential rotation and meridional flow are most reli- 
able at low to mid latitudes, this is where we focus our 
attention. These indicate a single circulation cell near 
the surface up to latitudes of at least 50-60°. According 
to the gyroscopic pumping equation (0, the presence 
of one or more high-latitude counter-cells may signify a 
change in the sense and possibly the amplitude of the 
Reynolds stress, which may in turn signify a change in 
V c . However, giv en the small latit udinal variation of the 
solar irradiance (|Rast et al.ll2008|) . it is more likely that 
a change in T at high latitudes would be associated with 
a change in the efficiency factor e, the transport scale <5, 
the angular momentum gradient |V£|, or other contri- 
butions to T such as Maxwell stresses. These issues lie 
outside the scope of this paper. 

3.2. Estimates Based on Helioseismic Inversions 

In Sj3J]we gave a conservative estimate of | V£| by con- 
sidering only the uniform rotation component. A more 
precise estimate follows straightforwardly from helioseis- 
mic rotational inversions and is shown in Fig. [5J At 
mid-latitudes its magnitude is roughly 3 km s _1 , slightly 
larger than the value associated with the uniform rota- 
tion component, |V£o| ~ 2Vh ~ 2.7 km s^ 1 . This is a 
consequence of the sense of the differential rotation, pro- 
grade at the equator and retrograde at the poles, which 
enhances the cylindrically outward V£. Thus, the sim- 
ple lower limit expressed in (0 is still valid for the Sun 
when the differential rotation is taken into account. The 
same is true for any star with a solar-like differential ro- 
tation (equatorward VQ). 

An estimate for V m is harder to come by but we do 
have two pieces of reliable information from solar obser- 
vations that we will exploit. First, the mean meridional 
flow at mid-latitudes is poleward on average near the sur- 
face for r > 0.95i?, with an amplitude of about 10-20 m 
s _1 (31]). Second, solar structure models provide a reli- 
able measure of the density throughout the convection 
zone that is verified by h elioseismic structure inversion s 
to within a few percent ([Christensen-Dalsgaardl 12002). 
We now proceed to use these two foundations together 
with mass conservation to obtain a lower bound on V m . 

We begin by estimating the poleward mass flux near 
the surface at mid-latitudes 

M = 2tt / (pvg) rdr = 2ir V g / prdr , (13) 

J r s J r s 

where r s = 0.95R and we will take p from Model S (see 
£|3.1[) . In order to obtain a lower limit on V m , and thus, 
the convective velocity V c , we take the mass- weighted 
poleward flow near the surface to be V$ ~ 10 m s" 1 , 
which gives M~4.1x 10 21 g s -1 . 
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Fig. 3. — Magnitude of the deep meridional flow V m as a function 
of the turnaround radius ro/R estimated from Eq. (1141) for rj, = 
0.69-R (grey shaded area) and r;, = 0.71 (hatched area). Each area 
is bounded by the lower curve, given by M = 4.1 X 10 21 g s _1 and 
an upper curve where M is given by Eq. H13H with r s replaced by 
ro and Vg = iOms" 1 . 

We treat M as a known quantity and estimate the 
average (density-weighted) deep return flow V m based 
on mass conservation as 



T,„ ~- M { 2k I prdr 



(14) 



Here rt, is the radius below which the meridional flow 
is negligible and ro is the turnaround radius where the 
mean flow shifts from poleward to equatorward (see Fig. 
Q]and the associated discussion in £|3.1[) . Note that this 
does not preclude multiple cells in radius; regardless of 
the meridional flow profile, there must be a net equator- 
ward flow below r s to balance the poleward flow above r s . 
Note also that some component of this return flow must 
necessarily cross C contours (which are approximately 
cylindrical: see MH11), requiring a Reynolds stress to 
maintain it as expressed in Eq. ([2]). 

Figure [3] shows estimates of V m based on two different 
assumptions for the base of the main circulation cell, rf,. 
A robust lower limit for V m is obtained by setting M 
to be the value known for the near-surface flow above 
0.95.R. An upper limit follows if we assume that the 10 
m s _1 poleward flow persists for the entire region above 
the turnaround radius r > rg. Then the deep flow V m 
must be correspondingly stronger to balance the greater 
poleward mass flux. The actual meridional flow speed is 
likely to lie between these limits, as indicated in Fig. [3] 

We are now able to obtain quantitative limits on V c 
based on the expressions derived in jj3.ll Taking | V£| ~ 
3 km s _1 at mid latitudes from Fig. [2] and 5 from Eq. ([9} 
and substitute them into Eq. JJJ gives 



V c > 27 m s- 1 N~ 1/2 



Vrr, 



1 m s 



1/2 



(15) 



Thus, if the latitudinal flow speed V m lies somewhere 
between 0.1-10 m s _1 as suggested by Figure [3l then the 
lower limit for V c in Eq. (|15|) lies somewhere between 8.6 
- 86 m s _1 . 



The higher end of the range, ~ 86 m s , will apply 
near the surface (r > 0.95i?) where the meridional flow 
speed is known to be 10-15 m s _1 . The lower estimate of 
8.6 m s _1 only applies if the meridional flow speed is as 
slow as 0.1 m s" 1 near the base of the CZ, if it penetrates 
deeper than 0.69i?, or if the circulation profile is multi- 
celled (N c > 1). Lower values of V c are also possible if 
the meridional flow near the base of the convection zone 
is confined to a narrow layer in radius as described in 
gHQ Near the top of the CZ, Eqs. fllQj) and © give a 
more reliable but more conservative estimate of V c > 45 
m s _1 . 

Recall that these estimates are based on the latitudinal 
flow speed at mid latitudes. By contrast, the estimate 
given in Eq. (fT2|) is based on the radial flow at the equa- 
tor. Thus, for this estimate, we use the larger equatorial 
value of | V£| ~ 5 km s _1 indicated by Fig. [3J We also 
equate the upward mass flux at the equator with the 
poleward mass flux above r ~ 0.95i? inferred from local 
helioseismology. As noted above, this gives a value of 
M = 4.1 x 10 21 g s^ 1 . Substituting these values into Eq. 
1(121) yields 



V r > 30 m s- 1 



1 


p 


-1/2 


r 




0.008 g cm" 3 _ 




Lo.95i?J 



(16) 



As noted in i j3.1[ this estimate is most reliable at r ~ 
0.95i? above which the meridional flow is known to be 
poleward so we have a reliable estimate for M. The value 
of p = 0.008 g cm -3 used in Eq. (fTfJl) is based on Model 
S at r ~ 0.95i?. If we assume this outward mass flux 
extends down to the lower CZ, then Eq. (p~6|) implies 
V c > 9 m s" 1 at r - 0.75i? (where p - 0.14 g s" 1 ). 

Combining these seperate estimates, we conclude that 
the amplitude of convective velocities must be at least 30 
m s -1 in the upper CZ (r ~ 0.95i?) and at least 8 m s _1 
in the lower CZ (r ~ 0.75i?). These of course are lower 
limits so larger amplitudes are entirely possible and even 
likely. 

4. DISCUSSION 
4.1. Length Scales and Implications for Detectability 

The lower limits obtained for convective velocities in 
make no reference to the characteristic length scale 
at which the convective motions occur. In this section 
we address this issue and consider its implications for 
the detectability of deep convection by means of local 
helioseismology. 

Some insight into the characteristic length scale L c 
for deep convection can be obtained by considering the 
Rossby number R a — V C /(2QL C )- As is well known, R a 
quantifies the amplitude of advective momentum trans- 
port in a rotating reference frame relative to the Coriolis 
force; small values R a < 1 imply strong rotational influ- 
ence and large values weak. 

Numerical simulations of convection exhibit a pro- 
found change in the nature of convective angular mo- 
mentum transport as the Rossby number varies across 
unity. When the rotational influence is weak R a > 1, 
convective flows tend to conserve their angular mo- 
mentum locally, producing anti-solar differential rota- 
tion profiles in which the angular velocity increases to- 
ward the rotation axis, implying dil/d8 < in the NH 
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20111: Ecathcrs tone et al.|[2012f ). As the rotational influ- 



ence becomes stronger, the preferred convective modes 
tend to align with the rotation axis, and the resulting 
velocity correlations induce an equatorward angular mo 
mentu m transport by means of the convec tive Reynolds 
stress (|Bussell2002HMiesch fc Toomrd [2009). This gener- 
ally promotes solar-like angular velocity profiles in which 
the angular velocity gradient is equatorward (9f2 /d6 > 
in the northern hemisphere). It is currently an open ques- 
tion precisely where this transition lies but simulations 
suggest that it occurs for values of R„ somewhat less 
than unity, perhaps around 0.6 ±0.3 (Fe atherstone et all 
[Ml) . 

Although this insight is based on convection simula- 
tions (which, like any model, have limitations), there is 
good reason to suspect that it still applies in the ex- 
treme parameter regime of the solar interior. Consider 
a radial downflow convective plume that is accelerated 
in the upper thermal boundary layer of the Sun, near 
the photosphere. The Coriolis force operating on this 
plume will tend to deflect it in a prograde direction with 
a characteristic time scale of r ~ (2J7o) _1 and a radius 
of curvature r c — V c /(2f2o). If the vertical coherence of 
the plume in the absence of rotation is L c then we obtain 
r c/L c = R . 

To appreciate the significance of this, consider the 
equatorial plane (sin6> = 1) and set L c to be the depth 
of the convection zone. If R >• 1, the plume is deflected 
only slightly in a prograde direction before it reaches the 
base of the convection zone. This will induce a nega- 
tive (v' r v'^\ correlation, producing cylindrically inward 

angular momentum transport and an anti-solar differen- 
tial rotation profile. However, if r c is less than the depth 
of the CZ, then our idealized, ballistic plume will never 
make it to the bottom. The vertical dominance of the 
flow will be broken and the nature of the Reynolds stress 
will be profoundly altered. This argument can be readily 
generalized to any intrinsic vertical coherence length L c 
and any latitude, with an effective radius of curvature of 
r c j sin#. 

More generally, it is reasonable to argue that only mo- 
tions with R < 1 will possess a rotational influence 
strong enough to establish the solar differential rota- 
tion (for r < 0.95i?). The corresponding length scale 
L c ~ V c /(2H, Ro) can be expressed in terms of the spher- 
ical harmonic degree £ ~ 27rr/i c , as follows; 



t~^Ro<750 R 

v c 



30 m s 



(17) 



For the numerical estimate we have used Oq = 2.7 x 10 -6 
and r = 0.95R. Note that £ = 750 corresponds to a 
physical length scale of 5.5 Mm. 

Thus, setting R Q < 1 in Eq. ([T7]) suggests that the 
convective motions responsible for maintaining the solar 
differential rotation and meridional circulation must oc- 
cur at spherical harmonic degree no greater than 750, 
implying a length scale no less than 5.5 Mm. 

Note that the values used to obtain the numerical esti- 
mate in Eq. fTT|. namely r = OMR and V c > 30 m s" 1 
correspond to the base of the near-surface shear layer 



where the transition from large to small Rossby number 
apparently occurs (MH11). Note also that these values 
are almost certainly a conservative estimate. If the tran- 
sitional Rossby number is closer to 0.5, the convective 
velocity scale is closer to 50 m s _1 and the efficiency fac- 
tor in Eq. ([6|) is closer to 0.5, then the lower limit in 
Eq- p7|) assumes plausible values of £ ~ 160 and L c ~ 
25 Mm. One might also argue that convective motions 
below 0.95i? are likely to occupy scales no smaller than 
supergranules, which would imply £ < 130 and L c > 30 
Mm. 

These latter values are remarkably close to the the den- 
sity and pressure scale heights, H p and Hp, which pro- 
vide an independent estimate for L c by virtue of mixing 
length theory. According to Model S, H p ~ 20 Mm and 
Hp ~ 12 Mm at r = 0.95R, increasing steadily to 90 
and 57 Mm respectively at the CZ base. Interestingly, 
if we assume that the length scale of convective motions 
at r = 0.95i? is equal to the density scale height, L c ~ 
20 Mm, then this implies that the transition from strong 
to weak rotational influence that marks the base of the 
NSSL occurs at a Rossby number of about 0.28. 

More generally, mixing-length theory predicts that L c 
should increase with decreasing r, becoming larger near 
the base of the CZ. A naive application of Eq. (fl7|) sug- 
gests the opposite; inserting our estimate of V c > 8 m 
s _1 in the lower CZ would push the limit toward higher 
£. However, this is a misapplication of Eq. (|T7|) . Convec- 
tion simulations and mixing length theory suggest that 
the Rossby number should decrease toward the base of 
the CZ more rapidly than V c , implying a n increase in L r . 

As mentio ned in £[TJ recent work by lHanasoge et al.l 
(|2010l I2012D based on local helioseismic inversions sug- 
gests that the spectral amplitude of convective motions 
may be no more than 1 m s _1 on scales £ < 60 at a radius 
of r ~ 0.92-0.95 R. Smaller horizontal scales (higher £) 
lie beyond their detection limits at that depth. This is 
difficult to reconcile with our lower limit of 30 m s _1 but 
may be possible if deep solar convection occupies multi- 
ple disparate scales, with broad, weak upflows surround- 
ing very narrow downflows. This may well yield small 
spectral amplitudes for global-scale modes (£ < 60) while 
local velocities in downflows might be much higher. How- 
ever, this still poses significant challenges to our current 
paradigm for how solar mean flows are maintained. 

In summary, the lower limits on convective velocities 
obtained in fj3] based on the mai ntenence of mea n flows 
and the upper limits obtained bv lHanasoge et all (|2010l ) 
based on local helioseismic inversions are both consistent 
with the idea that the characteristic length scale of con- 
vective motions is comparable to the local scale heights 
Hp and Hp. This yields Rossby numbers less than unity 
in the deep convection zone where mean flows are es- 
tablished and Rossby numbers greater than unity in the 
near-surface shear layer, as suggested by solar observa- 
tions (MH11). 

Finally, we note that the estimated length scales dis- 
cussed here bode well for the future of global solar con- 
vection simulations. In order to adequately resolve a 
structure of size L c , a simulation should have a resolution 
of at least ~ 0.1L C . Thus, if J? G ~ 1, Eq. fP7|) suggests 
that a resolution as high as £ ~ 7500 (grid spacing 6 ~ 
0.65 Mm) could be required to capture the relevant dy- 
namics. This is beyond the current capabilities of global 
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models but recall that this is a conservative limit. For 
plausible values in the mid CZ of R D ~ 0.2, V c ~ 50 
m s- 1 and e - 0.5, Eq. (HTJ) implies I < 64 (L c > 65 
Mm). Thus, a resolution extending to I ~ 640 should 
be sufficient to capture the dominant physical scales. 
This is achievable now with current global simulations 
(jMiesch et al.|[200l . 

This is not to say that current simulations necessar- 
ily capture all of the relevant dynamics. But, to put it 
colloquially, they're beginning to approach the right ball- 
park. This suggests that we may be near a threshold in 
the sense that moderate increases in resolution, together 
with improved modeling of the surface boundary layer 
and the overshoot region, could yield significant advances 
in our understanding of solar convection and the mean 
flows it establishes. This may be achieved both through 
the availability of increasingly powerful computing re- 
sources and through improved numerical algorithms that 
achieve higher parallel efficiency. 

4.2. Consistency with Convection Models 

We emphasize again that the limits on the convective 
velocities deduced in M3. II and £)3.2I do not depend on 
any theoretical or numerical model other than the basic 
MHD equations and the dynamical balance expressed by 
Eq. ©. Thus, it is of interest to ask whether theoretical 
and numerical models of convection are indeed consistent 
with these limits. 

The short answer is yes; models of solar convection 
based on several disparate physical perspectives and 
modeling approaches are uniformly consistent with the 
ideas propoposed in this paper. This is demonstrated in 
Figure 01 

Shown in the Figure are results from global convec- 
tion simulations based on the ASH code (black line, see 
Weber et al. 2011 for more information about this par- 
ticular simulation), surface convection simulations based 
on the MURAM code (red line, see Vogler et al. 2005 
and Rempel et al. 2009 for a description of the model 
and Rempel 2011 for more details on this series of simu- 
lations), and a hybrid model that combines surface con- 
vection simulations (from the STAGGER code) with a 
deep extrapolation based on mixing-length theory (blue 
line, see Trampedach & Stein 2011 for more information). 
Note that both surface convection simulations presented 
here correspond to the quiet sun, with no active regions 
or flux emergence. 

Superposed on the various curves are horizontal lines 
representing the theoretical lower limit for V c obtained 
from Eq. (JT5J) with V m = 1.0, 10, and 25 m s" 1 . All 
curves lie above the first limit, implying that the con- 
vective motions are at least in principle strong enough to 
sustain a solar-like differential rotation with a meridional 
flow of order 1 m s _1 or less in the deep CZ. However, 
they are only strong enough to sustain a meridional flow 
speed > 25 m s _1 in the upper half of the CZ. So, the 
convective amplitudes are consistent with mean flows in- 
ferred from solar observations and they confirm the ex- 
pectation from mass conservation that the flow speed of 
the deep equatorward return flow is likely to be less than 
the poleward surface value of 10-20 m s _1 . Further- 
more, all curves satisfy the limits obtained from Eq. (|16|) 
of V c > 30 m s _1 and V c > 9 m s -1 in the upper and 
lower CZ respectively and all curves satisfy the limit of 
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Fig. 4. — Comparison of theoretical lower limits for the convec- 
tion amplitude (dotted lines) with numerical and theoretical mod- 
els of convection (solid lines) . The three dotted lines are obtained 
from Eq. (1 1 51 ) with V m = 1-0, 10 and 25 m s _1 as indicated. The 
black line is obtained fro m a global convect ion simulation with the 
ASH code (described in Weber et al. 2011) and the brown line is 
from a simulation of sur f ace convection done with the MURAM 
code d Vogler et al.l 120051 : [Rempel et "all 120091) . The blue line is 
from a composite model combining surface convection simulations 
with the STAGGER code (r > 0.97J?) with a mixing length model 
(r < 0.97R), the latter calibrated to give th e same entropy .jump 
as th e simulation and scaled for continuity (Trampedach Sz Stein 
120111 1. 

V c > 45 m s _1 for r > 0.95i? obtained from eqs. ([7|) and 

Further confirmation of the consistency of these mod- 
els comes from mean-field models of the s olar diffe r ential 
rotat ion and meridional circulation by iRempell (|2005l 
2006). Here the convective amplitudes are essentially 
prescribed a priori by means of the imposed Reynold 
stress, modeled as a turbulent diffusion plus a A-effect. 
This is in contrast to the convection simulations shown 
in Figure [4] where V c is a product of the simulation. 
However, these mean-field models can be calibrated to 
produce solar-like mean flows so estimates of V c can be 
obtained by selecting the optimal transport coefficients. 
Here we are interested in particular in the A-effect, which 
is the non-diffusive component of the Reynolds stress ten- 
sor responsible for maintaining the differential rotation 
and, ultimately, the meridional circulation by means of 
gyroscopic pumping f i)2.1[) . 

In Rempel's models, the amplitude of the A-effect is 
given by Ao^t^o where Ao is a non-dimensional coeffi- 
cient of order unity and v t is the turbulent viscosity. Re- 
lating this to the convective velocity as in Eq. ([6]) yields 
V c ^ (Aot^Slo) 1 / 2 for e < 1. Solar-like mean flows are 
generally obtained with Ao ~ 0.8 and v t ~ 3-5 xlO 12 
cm 2 s _1 in the bulk of the CZ. Near the base of the CZ, 
v t drops sharply by more than an order of magnitude. 

The value of v t ~ 3-5 xlO 12 cm 2 s _1 in the upper CZ 
implies V c > 25-35 m s . This is comparable to our 
estimate of V c > 30 m s _1 based on Eq. (fT2|) . However, 
this lower limit lies below the limit of V c > 86 m s _1 in 
the upper CZ based on Eq. (fTBf and V c > 45 m s~ x based 
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on Eqs. and (jTUJ) . This implies that the effective value 
of (5 in the upper CZ is somewhat smaller than that given 
by Eqs. © and (TTU)) . Since v t is nearly constant in the 
upper CZ, this smaller value of 5 can be attributed to the 
density gradient, which factors into the Reynolds stress 
as indicated in Eq. Similarly, a strong convergence 
of the angular momentum flux (5 <C 1) near the base of 
the CZ (associated with the sharp drop in v t ) sustains 
an equatorward flow of a few m s _1 despite the low value 
of v t . 

The value of v t used by Rempel was deliberately chosen 
to be relatively small in order to make it more compat- 
able with the small value of the turbulent magnetic dif- 
fusivity rjt ~ 10 11 cm 2 s -1 in the mid CZ, needed for the 
operation of the advection-dominated flux-transport dy- 
namo (34.31) . In many mean-fi eld theories such as that of 
iKitchatinov fc Rtidigerl (|200l . v t is significantly larger, 
often exceeding 10 13 cm 2 s _1 as predicted by mixing- 
length theory. 

The consistency among these disparate models is re- 
markable. The surface convection simulations (blue and 
red curves in Fig. [4| make good contact with photo- 
spheric observations of convective amplitudes but they 
do not address mean flows; differential rotation and 
meridional circulation lie outside the scope of the mod- 
els. Meanwhile, the global convection simulations self- 
consistently produce a solar-like differential rotation and 
the convective velocity amplitudes roughly match the 
surface convection simulations (to within about 20%) in 
their small region of overlap near r ~ 0.97-R. Further- 
more, the global (ASH) simulation satisfies the constraint 
in Eq. ([7]) with i5 ~ 0.25 and meridional flow speeds are 
consistent with Figure EH with values of a few m s _1 
near the base of the CZ and ~ 1 0-20 m s " 1 near the 
surface. The mean- field models bv lRempell ((20051 [20061 
produce solar-like mean flows with Reynolds stress am- 
plitudes somewhat lower than suggested by the convec- 
tive models but they are consistent with an efficiency fac- 
tor of e < 1. Furthermore, they are still consistent with 
the velocity estimates put forth in Sj3]when one takes into 
account the relatively sharp gradients in the imposed an- 
gular momentum flux (small S), particularly at the base 
of the CZ. 

We also note that our lower limit for V c of 8 m s _1 
at the base of the CZ is consis tent with the upper limit 
of V c < 50 m s" 1 obtained by llsik fc Holzwarthl ([2009T ) 
based on numerical models of the interaction between 
thin flux tubes and convective flows. Larger convection 
amplitudes disrupt the flux storage, promoting buoyancy 
instabilities on time scales shorter than the dynamo am- 
plification time of a few years. 

4.3. Implications for Flux- Transport Dynamo Models 

Arguably one of the most successful (and certainly one 
of the most popular) current paradigms for modeling 
the origin of the solar activity cycle is the flux-transport 
dynam o. For recent re vi ews see iDikpati fe Gilmanl 
(2009) an d iCharbonneaul ([2010D and for further de - 
tails see IWang et al.l (11991): iChoudhuri et all (119951): 
IDikpati fe Charbonneaul (H999I ): IDikpati fe Gilmanl 
(1200 1L I2006D: iKiiker et"aT1 (120011) : iNandv fe Choudhuril 
120011): iBonanno et ail 120021): iRempell 12 006): 
Uouve fe Brunl (|2007T ): I Jiang et all (j2007l) : lYeates et al.l 



: IGuerrero et al.1 d2009l): lM unoz-J aramillo et al.l 
1201 ID : iHotta fc Yokovamal (|2010l ). and many 



Flux-transport (FT) dynamo models are mean-field 
models that solve the axisymmetric MHD induction 
equation, typically in the kinematic limit. Like other 
mean-field models, they involve specified prescriptions 
for magnetic field generation and transport by non- 
axisymmetric motions such as convection and flux emer- 
gence that lie outside the scope of the model. Their 
defining characteristic that sets them apart from other 
mean-field models is that the mean meridional circula- 
tion plays an essential role in transporting magnetic flux 
and in thereby regulating the cycle period. 

Of particular importance is the direction and speed 
of the meridional flow near the base of the convection 
zone. There, an equatorward flow with a speed of 2-4 
m s _1 provides a robust mechanism for producing solar- 
like butterfly diagrams, whereby the mean toroidal field 
(taken as a proxy for sunspots) migrates toward the equa- 
tor on a time scale of about 11 years. 

Most FT models are also Babcock-Leighton models in 
which the emergence and subsequent dispersal of active 
region flux in the photosphere acts as a source for mean 
poloidal field (e.g. lCharb onncau 2010). If this is the prin- 
cipal source of poloidal flux and if the mean toroidal 
field is generated near the base of the convection zone 
as in most models, then the dynamo requires a trans- 
port mechanism in order to operate. The two princi- 
pal transport mechanisms considered in the literature are 
the meridional circulation and turbulent transport, the 
latter typically represented as a turbulent diffusion or 
magnetic pumping. Thus, Babcock-Leighton FT mod- 
els may be further classified as advection-dominated or 
diffusion-dominated depending on w hich of these trans- 
port mechanisms plays a larger role (jYeates et al.l 120081 : 
IDikpati fc GilmarJl2T)09l ). 

In short, the meridional circulation in advection- 
dominated FT dynamo models serves two roles. First, 
it regulates the cycle period through the equatorward 
transport of toroidal flux near the base of the CZ and 
second, it couples the spatially separated source regions 
for mean poloidal and toroidal flux. These two roles are 
distinct but they are not independent. For example, if 
the transport mechanism coupling poloidal and toroidal 
sources is too efficient, it can "short-circuit" the dynamo, 
reducing the cycle period by decreasing the effective path 
length at the base of the CZ over which the equatorward 
advection of toroidal flux operates. 

The advection-dominated regime generally works well 
in the sense that it compares well with solar observa- 
tions. Imposing a single-celled meridional flow in ra- 
dius (N c — 1) with a poleward flow speed of 10-20 m 
s _1 near the surface as indicated by observations (SJTJ) 
and an equatorward return flow of 2-4 m s^" 1 near the 
base of the convection zone (consistent with Fig. [3]) gen- 
erally produces solar-like magnetic cycles with a dura- 
tion of about 11 years when used in conjuction with 
a solar-like differential rotation profile and a Babcock- 
Leighton source. More subtle issues such as dynamo 
parity, saturation, and cycle modulation have generally 
been handled through minor variations on this basic 
paradigm. For example, dipolar parity can be promoted 
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by an additional source of p oloidal flux deeper in the 
CZ (pikpati fc Gilmanl 120011 ) or by enhanced turbulen t 
diffusion in the surface layers ()Hotta &: Yo kovama 20101 ) . 

Although their empirical success and simplicity is 
rather compelling, the principle problem with FT dy- 
namo models has always been the theoretical justifica- 
tion of the advection-dominated regime. In order for 
the meridional circulation to dominate the flux trans- 
port, convective transport must be relatively inefficient. 
Advection-dominated FT models typically model the lat- 
ter as a turbulent diffusion and set the amplitude of 
the turbulent diffusivity rjt in the mid convection zone 



to be 



10 1 



cm s 



or less. This is at least two or- 



ders of magnitude smaller than estimates of r\t based 
on mixing- length theory and convection simulations; see 
iMunoz-Jaramillo et al.l (|2011f) for a more detailed discus- 
sion of the problem. 

Our estimates of the convective velocity and length 
scales in <J3] provide a measure of the convective trans- 
port that is independent of mixing length theory and 
convection simulations. Moreover, these estimates are 
directly linked to the amplitude and profile of the dif- 
ferential rotation and meridional circulation, which are 
essential ingredients in all FT dynamo models. Thus, to 
explore the implications of these estimates, we can take 
the actual mean flow profiles used in FT models as a 
starting point and then ask how strong the convective 
motions must be in order to maintain these flows. 

The essence of the problem can be appreciated merely 
from equation (J8|). This suggests that the convective ve- 
locity must be of order the geometric mean between the 
rotational velocity and the meridional flow velocity, pro- 
vided that the transport scale 5 is not inordinately small. 
As noted there, this implies V c > 47-106 m s~ x for the 
meridional flow speeds of 2-10 m s _1 typically used in 
FT dynamo models. If we then assume that the size scale 
of the convection L c is of order the density scale height of 
60 Mm in the mid CZ, we obtain r] t ~ V c L c /3 ~ 4 x 10 11 - 
2 x 10 12 cm 2 s _1 . This is significantly higher than the 
values used in the lower CZ for many FT dynamo models. 
In what follows we provide a more detailed exposition of 
this result and its implications. 

Thus, we begin with a standard FT dynamo model 
with a single-celled profile (N c = 1) and a solar-like dif- 
ferential rotation. If the turnaround radius r$ is near the 
middle of the CZ, then 8 is given by Eq. (j9j and a lower 
limit on V c is given by Eq. (JT5J) . Setting V m ~ 2 m s _1 
near the base of the CZ implies V c > 38 m s _1 . In the 
upper CZ we can set V m ~ 10 m s _1 which implies V c > 
86 m s _1 . 

By combining these estimates of the velocity scale with 
the limits on the length scale L c discussed in §4.11 we can 
obtain an estimate of the turbulent magnetic diffusivity 
rjt ~ V c L c /3. Taking the very conservative limit of L c > 
5.5 Mm gives i] t > 7 x 10 11 cm 2 s _1 near the bottom of 
the CZ and r) t > 10 12 cm 2 s _1 near the top. Although 
this is somewhat smaller t han estimates based on mixing - 
length theory (~ 10 13 , see lMunoz-Jaramillo et af1l201 If ) . 
it is still nearly an order of magnitude larger than the val- 
ues typically used in advection-dominated FT dynamo 
models. Using an estimate for L c based on density or 
pressure scale heights increases rjt by nearly another or- 
der of magnitude ( 34.1[) . 



As noted in §3.11 these limits can be avoided if there is 
a strong convergence of the angular momentum flux near 
the base of the CZ which would effectively confine the 
equatorward return flow to a thin layer (this applies to 
the return flow of several m s _1 required for the operation 
of the FT dynamo; weaker equatorward flows may exist 
outside of this layer). In particular, the lower limit on 

V c in Eq. (|7|) scales as L t where L t is the width of this 
layer. Reducing V c enough to yield an rjt of 10 11 cm 2 
s _1 would require that the equatorward return flow be 
confined to a layer no wider than 



r m V c z 



1.7 Mm - 2 x 10" 3 i? 



(18) 



Here we have used Eq. (J7]) with V c — 3r]t/L c , r\ t = 10 11 
cm 2 s" 1 , L c = 5.5 Mm, V m -2m s" 1 , |V£| - 3 km 
s _1 , S — \f2L t /r rn and r m — 0.7R. Larger values of 
L c or V m would imply even stronger gradients (smaller 
L t ). Furthermore, since scales as e _1 [see Eq. 01a 
value of e — 0.2 as suggested by convection models f £|3.ip 
would reduce the estimate of L t in Eq. (fT8f by a factor of 
five, to 340 km. No current FT models employ such an 
extremely asymmetric meridional circulation profile and 
it is questionable whether efficient equatorward transport 
of strong toroidal flux concentrations could even occur in 
such a thin layer. 

Even if such a thin layer were to exist near the base 
of the CZ, V c and r\ would still have to be large enough 
in the upper CZ to satisfy Eqs. (fT0|) [with Eq. 0] and 
(fl"6)) . These suggest that V c must be at least 30 m s _1 
at r ~ 0.95, independent of the deeper structure and 
amplitude of the meridional flow. If L c > 5.5 Mm then 
this suggests r\ t should be at least 5.5 x 10 11 cm 2 s -1 in 
the upper CZ for any model with solar-like mean flows. 
This is in fact satisfied by many current FT models, but 
again, this is a conservative estimate. If L c > 30 Mm 
and e — 0.2, this limit becomes much more stringent at 
rjt > 6.7 x 10 12 cm 2 s _1 . In the lower CZ, for N c = I, 
Eq. ((Til) suggests rjt > 1.6 x 10 11 cm 2 s" 1 for L c > 5.5 
Mm and e ~ 1 or rj t > 2 x 10 12 cm 2 s" 1 for L c > 30 Mm 
and e — 0.2. Multiple cells in radius (N c > 1) could in 
principle help mitigate these limits on V Cl rj, and L t but 
they are known t o have an adverse effe ct on the operation 
of FT dynamos (|Jouve fc Brunll200l . 

Thus, we conclude that current advection-dominated 
FT dynamo models are not strictly self-consistent. In 
particular, the amplitude of turbulent transport they 
typically assume is not commensurate with the mean 
flows they employ. Again, this has been argued before 
based on mixing-length theory and convection simula- 
tions but here we demonstrate more generally that it is 
a direct consequence of Eq. {2} . 

A potential way out of this dilemma is by moving 
to the diffusion-dominated regime. This would entail 
using values of rj t > 10 12 cm s _1 throughout most 
of the convection zone. Such models have indeed had 
some success in modeling the solar cycle, and may 
actually do better in certain respects than advection- 
dominated models. Examples include reproducing the 
observed cycle amplitude-period relationship and trigger- 
ing g rand minima through variations in the meridiona l 
flow (|Jiang et al.ll2007t lYeates et alj[200l : lKaraldl2010h . 
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However, it can be a challenge for diffusion-dominated 
FT models (or any model in which the transport by tur- 
bulent diffusion or magnetic pumping is very efficient) to 
achieve cycle periods as long as 11 years. As noted above, 
efficient turbulent transport in a Babcock-Leighton dy- 
namo tends to short-circuit the region over which the 
dynamo operates and to thus decrease the cycle period. 
This is sometimes avoided by allowing the meridional cir- 
culation to extend well below the convection zone. How- 
ever, if the turbulent diffusivity there is low, then these 
models may suffer from the same problem described here 
(namely that the local amplitude of rj t is incommensurate 
with the local amplitude of the meridional flow). Alter- 
natively, one could place the turnover radius tq higher 
up in the CZ to achieve a slower meridional flow in the 
lower CZ and thus a longer period (see Fig. 3). Still, the 
downward turbulent transport cannot be so efficient that 
it suppresses the poleward migration of residual flux from 
active regions that is seen in photospheric observations. 

Alternatively, it has been argued that the quench- 
ing of turbulent transport by Lorentz force feed- 
backs m ay be able to sa l vage t he advection-dominated 
regim e (jGuerrero et al.l 120091 : IMunoz-Jaramillo et al.l 
1201 If ). However, it remains to be seen whether or not this 
is viable. Quenching would be most effective for strong 
toroidal fields near the base of the CZ as opposed to the 
relatively weak poloidal fields whose turbulent transport 
across the CZ could potentially short-circuit the advec- 
tion by the meridional flow. In other words, quench- 
ing is more likely to regulate the first role noted above 
for the meridional circulation in FT dynamos (equator- 
ward advection of toroidal flux) as opposed to the second 
(coupling of poloidal and toroidal sources), and it is the 
second which generally defines the advection-dominated 
regime. 

Furthermore, the arguments presented here demon- 
strate that a significant reduction in V c by Lorentz force 
feedbacks as captured by quenching mechanisms would 
have substantial consequences for the mean flows. For 
example, if we treat Vq, as fixed, Eq. ([8]) indicates that a 
reduction in V c by an order of magnitude may be accom- 
panied by a reduction of the meridional flow speed V m 
by two orders of magnitude (unless the efficiency factor 
e increases) . This follows from the concept of gyroscopic 
pumping discussed in £12.11 provided that the nature of 
the Lorentz forces is solely to suppress convection. How- 
ever, if the Lorentz force exerts its own mean torques, 
which is likely, then this should be included in the right- 
hand-side of Eq. ([2]), with corresponding modifications 
to Eqs. (O and ([8]). The time dependence should also 
be taken into account, potentially mitigating limits on 
7/ f for a given meridional flow speed. Even so, a self- 
consistent treatment of diffusivity quenching that takes 
into account its dynamical effect on the maintenance of 
mean flows could dramatically alter the operation of a 
FT dynamo by substantially modifying the meridional 
circulation and, to a lesser extent, the differential rota- 
tion. 

A related possibility is that turbulent transport at very 
high magnetic Reynolds numbers may simply be less effi- 
cient than suggested by crude, kinematic representations 
based on the concept of turbulent diffusion. Nonlinear 
processes such as the dynamical alignment of fields and 
flows that contribute to dynamo saturation may also sup- 



press turbulent transport in a way that is more subtle 
than a quenched local diffusion coefficient. This might 
require more sophisticated mean- field and/or MHD con- 
vection models to properly capture. Still another pos- 
sibility, of course, is that the solar dynamo does not 
follow the canonical Babcock-Leighton/Flux- Transport 
paradigm. 

5. SUMMARY 

In this paper we estimate the amplitude and scale 
of the convective motions responsible for maintaining 
the solar differential rotation and meridional circulation. 
This estimate is based only on the observed properties of 
the mean flows in the Sun together with three fundamen- 
tal physical premises grounded in the MHD equations. 

The first and most important of these three premises 
is represented by Eq. @. This tells us that, in or- 
der to achieve a statistically steady state, the angular 
momentum transport by the convective Reynolds stress 
must balance the advection of angular momentum by the 
meridional flow. This provides a direct link between the 
amplitude of mean flows and the amplitude of convective 
motions such that observations of the former can set con- 
straints on the latter. Note that this dynamical balance 
is to be understood in a time-averaged sense, filtering out 
solar cycle variations of meridional and zonal flows. 

The second physical premise we rely on is that of mass 
conservation [Eq. (p} in the Appendix]. This allows us 
to estimate the net equatorward flow speed in the deep 
CZ based on observations of the poleward flow near the 
surface as shown in Fig. [3l This result in turn allows us 
to extend the diagnostic power of our first premise to the 
deep convection zone where mean flows are established. 

The third premise is that of thermal wind balance, Eq. 
([5]). Although this is not used directly in our estimates 
for the convective velocity amplitude V c , it helps to jus- 
tify the concept of gyroscopic pumping discussed in ij2.ll 
In particular, when coupled with helioseismic rotational 
inversions, it suggests that the principle mechanism re- 
sponsible for maintaining the solar meridional circulation 
is the inertia of the differential rotation, as represented 
by the Coriolis force ( ^T^, . 

Together with these three physical premises, we use 
three observational foundations. The first is the solar 
differential rotation inferred from helioseismology. The 
second is a persistent poleward meridional flow from low 
to high latitudes for r > 0.95-R with an amplitude of 
roughly 10~15ms _1 . The third observational foundation 
is the mean density profile p(r) in the CZ. Although this 
is obtained from a solar structure model (Model S), we 
regard it as an observation because it is verified to within 
a few percent by helioseismic structure inversions. 

The results suggest that the amplitude of convective 
velocities in the upper CZ (r ~ 0.95i?) must be at least 
30 m s _1 in order to sustain the observed mean flows. 
Analogous limits in the lower CZ are less reliable due to 
uncertainties about the meridional flow but reasonable 
inferences suggest that convective amplitudes can be no 
less than about 8 m s _1 at r ~ 0.75. 

The existence of the near-surface shear layer (NSSL) 
provides a smoking gun that we can exploit to link 
these convective amplitudes to a size scale. In partic- 
ular, it suggests that the characteristic Rossby number 
crosses unity somewhere near the base of the NSSL at 
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r ~ 0.95i? (MH11). Together with the velocity limits 
obtained in this implies deeper convective motions 
can be no smaller than 5.5 Mm. This is of course a 
conservative limit. Surface convection simulations sug- 
gest that the size scale of convective motions progres- 
sively increases with depth, from the ~ 1 Mm granu- 
lation scales at the surface to ~ 30 Mm supergranula- 
tion scales by r ~ 0.97-R, due mainly to an increase 
in the density and pressure scal e heights (e.g. iRempell 
1201 U iTrampedach fc Steinll201lD . Still, the limit given 
here provides an independent constraint on the size of 
deep convection that depends only on the three physi- 
cal premises and observational foundations noted above 
(along with a fourth premise that R Q < 1 marks the base 
oftheNSSL). 

The lower limits on convective velocity and size scales 
reported here are consistent with the upper limits in- 
ferred from local helioseismology by Hanasog e et al.l 
(2010) provided that the characteristic size of the convec- 
tive motions at r ~ 0.95i? lies somewhere in the range be- 
tween 5.5-83 Mm. This corresponds to a range in spher- 
ical harmonic degree of 50 < £ < 750. This range may 
be too wide to provide a strong constraint on convection 
models but it can be narrowed somewhat if one requires 
that motions at and below r ~ 0.95i? be no smaller than 
supergranules, so L c > 30 Mm and i < 130. 

Since our work is concerned specifically with the rela- 
tionship between convective transport and mean flows, it 
has important implications for Flux Transport dynamo 
models f ^4.3[) . In particular, it suggests that advection- 
dominated FT models are not self-consistent in the sense 
that the assumed magnitude of convective transport is 
generally too low to account for the mean flows they re- 
quire to operate. This has been argued before based on 
estimates of V c and L c obtained from mixing length the- 
ory and convection models (cf. Fig 3). However, here we 
demonstrate that it is a more general consequence of the 
need to sustain the mean flows against their own inertia. 

The essence of the problem can be appreciated simply 
from Eq. ©. If (2S) 1 / 2 ~ 1 and V n > V m , then V c 
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must be significantly larger than V m . Thus, one would 
expect convection to dominate transport over the merid- 
ional circulation unless its spatial and temporal correla- 
tion scales are small. However, these correlation scales 
cannot be too small because if they were, the influence of 
rotation on the convection would not be strong enough to 
establish a solar-like differential rotation (i.e. R would 
be greater than unity; see £|4.1[) . 

This result casts some doubt on the advection- 
dominated flux-transport paradigm as a viable model of 
the solar cycle but it does not necessarily rule it out. An 
alternative possibility is that convective flux transport 
in the extreme parameter regimes of the solar interior is 
much less efficient than suggested by turbulent diffusion 
and that we have much to learn about reliably represent- 
ing this transport in mean-field dynamo models. 

The estimates reported here are independent of 
mixing-length theory and convection simulations but 
they give compatible valu es for V c and L c , suggesting 
internal consistency f £|4.2p . Furthermore, they suggest 
that the velocity and length scales responsible for main- 
taining solar mean flows cannot be drastically different 
from those currently exhibited by global convection sim- 
ulations. This bodes well for the future; as global con- 
vection models continue to move toward higher resolu- 
tion and improved representations of the upper and lower 
boundary layers, they should be able to capture the es- 
sential physics underlying the solar differential rotation 
and meridional circulation with increasing fidelity. 
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APPENDIX 

Baroclinicity as a Source of Differential Rotation 

In ^2.2\ we argued that baroclinic forcing cannot account for the observed sense of the solar meridional flow (poleward 
near the surface). In this Appendix we demonstrate further that baroclinicity alone cannot account for the solar 
differential rotation as inferred from helioseismology . Similar issues have been studied for decades within the context 
of planetary and stellar atmospheres ()Eliassenlll951t [Read 19861 rTassoullll978l ). 

The importance of baroclinicity for shaping the solar differential rotation profile is undisputed. Theoretical models, 
mean-field models, and global convection simulations all suggest that baroclinic forcing i s necessary to account for the 
conical nature of mid-la t itude fl surfaces in the solar CZ inferred from helioseismolog y (Kitchatinov k Rudigcrj |1995t 
iRobinson k Charil200ll : iBrun k Toonirel l20ljl iRempell [2001 iMiesch et al1[200l iBalbus et all 12009( 1 . However, this 
result should not be misinterpreted to attribute the solar differential rot ation entirely to baroclin ic forcing. 

Such a clarification is particularly timely in light of the recent work by IBalbus k Schaanl ()2012h . They demonstrate 
that the centrifugal distortion of the base of the convection zone can produce thermal gradients along isobaric surfaces 
that can in turn induce differential rotation through baroclinic forcing. This background shear may then interact 
with convection to produce the observed mean flows. We do not dispute this argument but again we caution against 
a potential misinterpretation of their results. In particular, we argue that baroclinicity alone cannot induce a net 
equatorward angular velocity gradient (d£l/d8 > in the northern hemisphere, or NH) throughout the CZ as exhibited 
by helioseismic rotational inversions. 

Our argument begins with the time-dependent zonal vorticity equation (j3J). Again we neglect the Reynolds stress, 
Lorentz force, and viscous diffusion. We proceed to consider a fixed entropy gradient d (S) jdQ and we ask what 
differential rotation such baroclinic forcing will produce when subject to a given initial condition Cl(r, 9, t = 0) = fii(r, 6) 
and v m (i = 0) = 0. 

As noted in W1.2\ a poleward entropy gradient will induce a clockwise meridional circulation cell in the NH 
(d (S) Jd9 < 0, (aty) > 0). Although it was not necessary for the arguments presented in S}2l we now adopt the 
anelastic approximation so the mean mass flux is divergenceless: 

V- (pv ro ) = (1) 

where p(r) is the background density profile, averaged over latitude, longitude, and time. Thus, the conservation of 
mass requires that the generation of zonal vorticity (w^) be associated with flows both toward and away from the 
rotation axis. This is crucial in understanding the significance of the results that follow. 

In order to address the generation of rotational shear, we must consider the time-dependent version of the zonal 
momentum Eq. ((T|). Another crucial realization required to appreciate the results that follow is that there is no 
baroclinic component of the net axial torque J- . Recall that T involves only the Reynolds stress, the Lorentz force, 
and the viscous diffusion (explicit expressions are given in MH11). If baroclinic forcing is to induce a differential 
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rotation, it must do so by means of the meridional flow. Thus, we can set T = 0, which yields (MH11) 

dC 

p— = - (pv m ) • V£ . (2) 

In Eq. @ we have again used the anelastic approximation so we have replaced p with ~p. 

Two things are immediately apparent from Eq. First, for a given meridional flow (pv m ), this is identical to the 
equation for passive advection of a scalar field. Second, the only steady solution is one in which (pv m ) • V£ = 0. Such 
a steady state can be achieved in one of two ways. Either the meridional flow must vanish (v m ) = or the specific 
angular momentum C must be constant on streamlines. The latter case, L = A 2 il constant on streamlines, would 
necessarily be associated with an anti-solar differential rotation profile, such that the poles would spin faster than the 
equator (dil/09 < in the NH). We will return to the former case, that of (v m ) = below. 

Now return the time-dependent problem outlined above in which we follow the response of the mean flows to a 
specified baroclinic forcing d (S) / 89. Let us assume for simplicity that the initial angular velocity profile is cylindrical, 
so Qi — f2j(A). As noted below, this is not a necessary assumption but it serves well to illustrate the main point. 
Furthermore, we will assume that the initial angular momentum gradient is directed away from the rotation axis, so 
dCi/dX > 0, where Li = X 2 fli. This is true for the current Sun, it is true for a uniform rotation (fij = £Iq), and 
we suspe ct that it is tr ue in general for stars since the alternative, dCi /dX < is unstable according to the Rayleigh 
criterion (|TassouH[l978). Note that the case V£ = is a fixed point; if this were the initial state, meridional circulation 
would not influence the zonal flow and baroclinicity would not induce a differential rotation. 

Consider a volume V defined by r < R and A < Ao, where R is the solar surface and Ao is a fiducial cylindrical radius 
that can lie anywhere between zero and R. Integrating Eq. @ over V and using Eq. ([1} yields 

^ = £ L (pv m ) -dS = 2irX J + £ (pvx) dz (3) 

where 

Ly= / pCdV (4) 



JV 

is the total angular momentum within V and z± — ±(R 2 — A 2 ,) 1 / 2 . We have assumed that there is no flow across 
r = R. 

At the initial time t — the f2 profile is cylindrical so C can be pulled out of the integral on the right-hand-side 
(RHS). Mass conservation then implies that the RHS vanishes and the rate of change of Ly vanishes. However, the 
analogy with the passive scalar noted above makes it clear that this cannot hold indefinitely Our intuition tells us 
that if dCi I dX > then the advection of angular momentum into our volume must increase Ly over time. 

This apparent inconsistency can be resolved if we consider a Taylor expansion at early times 



L v {t) = L v (0) + t 



+ t 2 



d 2 L x 



dt 

Applying a time derivative to (j3]) and substituting in ([2]) yields 

d 2 L v 



dt 2 



(5) 



dt 2 



t=o 



2-tA / + (pvx) ^dz = 2ttA f + p (v x ) 2 ^ dz . (6) 



In obtaining this result we have neglected the time dependence of (v m ). Thus, the second derivative of Ly is positive 
definite at t = 0. Together with ([5]) this implies that the angular momentum within V will increase with time. 

Put another way, this implies that any persistent meridional circulation will tend to spin up the poles relative to the 
equator. This is true regardless of the (nonzero) amplitude or profile of the circulation and regardless of the source 
of the circulation, whether it be baroclinic in nature or due to some other meridional forcing. Furthermore, it holds 
for any arbitrary initial C profile that is stable in the sense that dCi/dX > 0. In this general case, C surfaces would 
still intersect R at two locations and V would be defined relative to the surface C = C$ such that r < R and C < Cq. 
Then Eq. (j6|) would still hold with dCi/dX replaced by | V£|, (v\) 2 by ((v m ) •n) 2 (where h is the unit vector normal 
to S) and the integration would proceed over the surface S. The early time dependence would still be given by ([5]) 
with (dL v /dt) t=0 = and (d 2 L v / dt 2 ) t=0 > 0. 

This argument can also be readily generalized to a spherical annulus. If there is a radius r& below which the 
meridional flow (v m ) becomes negligible, then one can define the volume V to be bounded from above and below by R 
and rb, extending poleward of an £ isosurface that extends from R to (confined either to the northern or southern 
hemisphere). Then the conclusion is unaltered; a baroclinically-driven meridional flow will increase the total angular 
momentum in the polar cap V over time, provided only that dCi/dX > 0. 

Spin-up of the poles can be avoided if the Rossby number is small, which is a good approximation for the Sun. Here 
the relevant Rossby number is that based on the meridional flow R" 1 = V m / (2Q,oR). For R™ -C 1, meridional flows are 
redirected into zonal flows essentially immediately (on the time scale of a rotation period), so the steady-state version 
of Eq. @ can be satisfied with (v m ) = 0. At the same time, a differential rotation will be quickly established that 
obeys thermal wind balance, Eq. ([5]). 
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This is indeed a viable way to generate an axial shear dQ/dz and may play an important role in the Sun. However, 
it cannot spin up the equator relative to the poles, as seen in the solar convection zone. In the limit R" 1 <C 1, there is 
no angular momentum transport across A surfaces. Rather, at a given cylindrical radius A, any deceleration of f2 for 
some range of z must be balanced by an acceleration at some other z. This follows directly from the conservation of 
mass and angular momentum as outlined above. 

This implies that a solar-like differential rotation characterized by an equatorward Vf2 in the CZ can only be 
produced by baroclinic forcing if there is a corresponding poleward Vf2 somewhere below the CZ. We cannot rule this 
out from helioseismic rotational inversions but there is currently no evidence for it. 

We emphasize that the thermal wind equation (J5J) possesses a geostrophic degeneracy in that a given entropy gradient 
d (S) /d9 is compatable with an infinite number of ft profiles, some solar-like {dVl/dO > in the NH) and some anti- 
solar. This arises because one can take a solution of this equation f2 s (r, 9) and add an arbiratry cylindrical component 
f2'(A) and the sum is also a solution. We argue that the convective Reynolds stress is necessary to break this degeneracy 
and establish a solar-like ft profile. 

In any case, it is clear that there must be some non-baroclinic forcing in the solar CZ in order to account for the 
observed mean flows. As noted in <JTJ solar observations indicate a persistent poleward circulation in the solar surface 
layers and mass conservation requires a net equatorward flow deeper down. Since C isosurfaces are nearly cylindrical 
throughout the CZ (MH11), this necessarily implies flow across C contours that must be balanced in a steady state 
with a nonzero J 7 , as expressed in Eq. ([1}. 



